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The four major asymptotic level density laws of random matrix theory may all be show¬ 
cased though their Jacobi parameter representation as having a bordered Toeplitz form. 
We compare and contrast these laws, completing and exploring their representations in 
one place. Inspired by the bordered Toeplitz form, we propose an algorithm for the finite 
moment problem by proposing a solution whose density has a bordered Toeplitz form. 

Keywords : Finite moment problem; Infinite random matrix theory; Jacobi parameters; 
Toeplitz matrix. 


1. Introduction 

Consider the “big” laws for asymptotic level densities for various random matrices: 

Wigner semicircle law 
Marclienko-Pastur law 
Kesten-McKay law 
Wachter law 

In raw form, these laws (Table [lj appear as somewhat complicated expressions 
involving square roots. This paper highlights a unifying principle that underlies 
these four laws, namely the laws may be encoded as Jacobi symmetric tridiagonal 
matrices that are Toeplitz with a length 1 boundary. 

This suggests that some of the nice properties of the big laws are connected 
to this property, and further suggests the importance of the larger family of laws 
encoded as Toeplitz with length k boundary, known as “nearly Toeplitz” matrices. 
This motivates the two parts of this paper: 

(1) We tabulate in one place key properties of the four laws, not all of which can be 
found in the literature. These sections are expository, with the exception of the 
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Measure 

Support 

Parameters 

Wigner semicircle 
yjA-x 2 

dfJ'WS = 2it ax 

Iws = [±2] 

N/A 

Marchenko-Pastur 
d-UMP ~ ^ (yX+ ~ X ^ “ A -) dx 

Imp = [A_, A + ] 

A± = (1 i^A) 2 , A > 1 

Kesten-McKay 

v^/A(v — 1) - x 2 

dfiKM = dx 

2tt(v z — x z ) 

Im = [±2i/v — 1] 

v > 2 

Wachter 

{a + b)^{n + -x){x-nJ) 
d <* w ~ 2»(1 - x) dX 

hv = [fj,-,n+] 

/ Vb± \J a{a + b — 1) A 

a + b ) ’ 

a, b > 1 


Table 1. Random matrix laws in raw form. The Kesten-McKay and Wachter laws are related by 
the linear transform (2a:yy ac j 1 ^ er — l)f = - T Kesten-McKay an< ^ a = b = v/2. 


as-of-yet unpublished Wachter moments, and the Kesten-McKay and Wachter 
law Jacobi parameters and free cumulants. 

(2) We describe a new algorithm to exploit the Toeplitz-with-length-fc boundary 
structure. In particular, we show how practical it is to approximate distribu¬ 
tions with incomplete information using distributions having nearly-Toeplitz 
encodings. 

Studies of nearly Toeplitz matrices in random matrix theory have been pioneered 
by Anshelevich m- 

Historically, the Wigner semicircle law is the most famous. The weight function 
is classical, and corresponds to Chebychev polynomials of the second kind. It is the 
equilibrium measure [1] for Hermite Polynomials and the asymptotic distribution 
for Gaussian or Hermite ensembles (GOE, GUE, GSE, etc.). None of the other 
weight functions are classical, but they are all equilibrium measures for classical 
polynomials. The second most famous law is the Marchenko-Pastur law. It is the 
equilibrium measure for Laguerre Polynomials and is the asymptotic distribution 
for Wishart matrices or Laguerre ensembles. The Kesten-McKay law, described in 
m , is the equilibrium measure for Gegenbauer Polynomials. It is not commonly 
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included among the Wigner semicircle, Marclienko-Pastur, and Wachter laws, but 
we believe that it merits inclusion on account of its place in the upper-right box 
in Table [2] The Wigner and arcsine distributions are special cases of the Kesten- 
McKay distribution. The Wachter Law generalizes to Jacobi Polynomials. They 
describe sections of random unitary matrices, MANOVA matrices, and general Ja¬ 
cobi ensembles. 

The laws have been encoded in many formats over the years. Wigner’s earliest 
work encoded the semicircle law through its exponential generating function, the 
Bessel function. For many decades the Stieltjes or Cauchy transforms (Table[3]) have 
been valuable. Free probability has proposed the R-transform and S-transform as 
efficient encodings for these laws (Table [5]). In this note, we list each representation 
of each distribution, but we focus on the Jacobi parameter representation. 

Other laws may be characterized as being asymptotically Toeplitz or numerically 
Toeplitz fairly quickly, such as the limiting histogram of the eigenvalues of {X/yJm-\- 
ftiy (X / s/m + pi), where X is m x n, n is 0 (m), and m —> 00 (Figure [3]). 

This nearly Toeplitz property inspires an underlying approximation concept. 
Instead of simply truncating a Jacobi matrix, we can construct an infinite Toeplitz 
matrix that is Toeplitz on all but an initial finite set of rows and columns. In 
turn, we apply this approximation idea to the moment problem in Section [5j given 
a finite set of moments, we can use the Lanczos iteration to find corresponding 
Jacobi parameters, put those parameters in a continued fraction, and take its inverse 
Stieltjes transform to find a smooth, compactly-supported distribution with the 
correct moments. Instead of using the moments to find the Jacobi parameters we 
can also use a discretization of the measure. In :lH], the Cauchy transform and 
its relationship to continued fractions of depth 2 are discussed, an idea that is 
generalized here in Table [9j 

Mathematical investigations into weight functions with Toeplitz or asymptoti¬ 
cally Toeplitz Jacobi parameters may be found in the work of two mathematicians 
with coincidentally similar names: Geronimus [81 and Geronimo 1, It seems 
known that the algorithm may have issues, perhaps reminiscent of the Gibbs phe¬ 
nomenon of Fourier analysis, in that atoms may emerge. 

In retrospect, a good part of this paper may be found explicitly or implicitly in 
the work of Anshelevich [2]. Nonetheless, as we began to form the various tables 
and noticed the nice Catalan and Narayana properties, and especially the little box 
in the upper right in Table [2j we realized that in addition to the algorithm, there 
were enough ideas that were below the surface that we worked out for ourselves 
and wished to share. 


2. The Jacobi Symmetric Tridiagonal Encoding of Probability 
Distributions 

All distributions have corresponding tridiagonal matrices of Jacobi parameters. 
They may be computed, for example, by the continuous Lanczos iteration, described 
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in .191 P-286] and reproduced in Table [8] 

We computed the Jacobi representations of the four laws providing the results 
in Table[2j The Jacobi parameters (a* and p\ for i = 0,1,2,...) are elements of an 
infinite Toeplitz tridiagonal representations bordered by the first row and column, 
which may have different values from the Toeplitz part of the matrix. 


a 0 

Po 

Po 

ai Pi 


Pi an Pi 


Pi OL\ Pi 


Pi ai 



c*o = ai 

a 0 ± a 1 

Po — Pi 

Wigner semicircle 

Marchenko-Pastur 

Po^Pi 

Kesten-McKay 

Wachter 


Measure 

an 

Ctn, (n > 1 ) 

Po 

Pm {n > 1 ) 

Wigner Semicircle 

0 

0 

1 

1 

Marchenko-Pastur 

A 

A + l 

\/A 

VX 

Kesten-McKay 

0 

0 


\/v — 1 

Wachter 

a 

a 2 — a + ab + b 

Vab 

\Jab(a + b — 1) 

a + b 

(a + b) 2 

0 a + 6) 3 / 2 

(a + b) 2 


Table 2. Jacobi parameter encodings for the big level density laws. Upper left: Symmetric Toeplitz 
Tridiagonal with 1-boundary , Upper Right: Laws Organized by Toeplitz Property, Below: Specific 
Parameter Values 


Anshelovich [2] provides a complete table of six distributions that have Toeplitz 
Jacobi structure. The first three of which are semicircle, Marchenko-Pastur, and 
Wachter. The other three distributions occupy the same box as Wachter in Ta¬ 
ble [2] Anshelovich casts the problem as the description of all distributions whose 
orthogonal polynomials have generating functions of the form 

OO 

J2 P n(x)z n = 

n—0 


1 — xu(z) + tv{z) ’ 
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which he calls Free Meixner distributions. 

He includes the one and two atom forms of the Marchenko-Pastur and Wachter 
laws which correspond in random matrix theory to the choices of tall-and-skinny 
vs. short-and-fat matrices in the SVD or CS decompositions, respectively. 

3. Infinite RMT Laws. 

This section compares the properties of all four major infinite random matrix theory 
laws, the Wigner semicircle law, the Marchenko-Pastur law, the Kesten-McKay law, 
and the Wachter law. 

We state the four laws of infinite dimensional random matrix theory and their 
intervals of support, pws on Iws, Bmp on Imp, Rkm on I km, and pw on Iw cor¬ 
respond to the Wigner semicircle, Marchenko-Pastur, Kesten-McKay, and Wachter 
laws. These laws are originally credited to m, M, Ha, and m- See Tabled 


Measure 

Cauchy Transform 

Wigner Semicircle 

1 

CN 

1 

Marchenko-Pastur 

1 — X + z ~ y/0- — A + z) 2 — 4 z 

2z 

Kesten-McKay 

(v — 2 )z — V\jA(l — v) + z 2 

2(v 2 — z 2 ) 



Wachter 

1 — a + (a + b — 2 )z — y/{a + 1 — (a + b)z) 2 — 4a(l — z) 

2z(l — z) 


Table 3. Cauchy transforms . 


We can also write down the moments for each measure in Table [4j for Wigner 
and Marchenko-Pastur see [5], for Kesten-McKay see m, and for Wachter see 
Theorem 6.1 in the Section 6. Remember the Catalan number C n = —( 2n ) and the 
Narayana polynomial N n (r) = J2j=i N n,jr 3 , where N n j = excepting 

N 0 (r) = 1. The coefficients of ib(l — u) 71 / 2- - 7 in the Kesten-McKay moments form 
the Catalan triangle. We discuss the pyramid created by the Wachter moments in 
Section [4j 
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Measure 

Moment n 

Wigner Semicircle 

C n /2 if n is even, 0 otherwise 

Marchenko-Pastur 

N n (X) 


Kesten-McKay 


n/2 

E 

3 =1 


n-j\ / j 


n/2 J \n — j 


V-* (y — l) n / 2 1 if n is even, 0 otherwise 


Wachter 


a 

a + b 


n —2 

- ( a + 6 ) E 

l=o 


2j+4 


y'aCa + fr- 1 ) \ 

a + b J J+1 ^0(0 + b — 1 ) 


Table 4. Moments 


Inverting the Cauchy transforms and subtracting 1/w, computes the R- 
transform, see Tabic [5] If there are multiple roots, we pick one with a series expan¬ 
sion with no pole at w = 0. 

The free cumulants K n for each measure appear in Table [6] by expanding the 
I?-transform above (the generating function for the Narayana polynomials is given 
by US], the generating function for the Catalan numbers is well known). 

It is widely known that the Catalan numbers are the moments of the semicircle 
law, but we have not seen any mention that the same numbers figure prominently 
as the free cumulants of the Kesten-McKay Law. The Narayana Polynomials are 
prominent as the moments of the Marchenko-Pastur Law, but they also figure 
clearly as the free cumulants of the Wachter Law. There are well known relation¬ 
ships, involving Catalan numbers, between the moments and free cumulants of any 
law m, but we do not know if the pattern is general enough to take the moments 
of one law, transform it somewhat, and have them show up in the free cumulants 
in another law. 

We compute an S-transform as S(z) = R~ 1 (z)/z. See Table[5] 

Each measure has a corresponding three-term recurrence for its orthonormal 
polynomial basis, with g_i(x) = 0, qo(x) = 1, /3_i = 0, and for n > 0, q n +i(x) = 
{(x — a n )q n (x) — P n -iq n -i(x))//3 n . In the case of the Wigner semicircle, Marchenko- 
Pastur, Kesten-McKay, and Wachter laws, the Jacobi parameters a n and j3 n are 
constant for n > 1 because they are all versions of the Meixner law [2] (a linear 
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Measure 

-R-transform 

R-transform 

Wigner Semicircle 

w 

1 

Marclrenko-Pastur 

A 

z — A 

1 — w 

2 2 

Kesten-McKay 

—V + Vy/l + 4 W 2 

V 

2 w 

V 2 — Z 2 




Wachter 

—a — b + w + \J (a + b) 2 + 2 (a — b)w + w 2 

a — az — bz 

2 w 

z 2 (z - 1) 


Table 5. /^-transforms and S'-transforms computed as S(z) = R 1 (z)/z . 


Measure 


Wigner Semicircle 

fin, 2 

Marchenko-Pastur 

A 

Kesten-McKay 

(—l)( n_2 )/ 2 uC '(„_ 2 )/2 if n is even, 0 otherwise 

Wachter 


v a) (a + b) 2n+1 


Table 6. Free cumulants. 
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transformation may be needed). The Wigner Semicircle case is given by simplifying 
the Meixner law in [Q, and the Marchenko-Pastur, Kesten-McKay, and Wachter 
cases are given by taking two iterations Lanczos algorithm symbolically to get a± 
and /3i. See Table [2] 


Measure 

q n {x), n> 1. 

Wigner Semicircle 

Ml) 

Marchenko-Pastur 

A<—>/>(x »U„n i^ 1 ) A"^£/ n _ 2 (-A=i) 

Kesten-McKay 

(v - 1)("— 1 V^xUn-! (^y) - v(v - l)("- 2 )/ 2 t/„_ 2 (^) 

Wachter 

( , a \ f \/ ab(a+b— 1 ) ^ TJ f —b—a(a+b—l)+(a+b) 2 x\ 

V a + b J { (“+ b ) 2 J n_1 V 2yjab(a+b-l) J 

a+b ( \/ ab(a+b— 1 ) \ „ / — b—a(a+b— l) + (a+b) 2 x \ 

a+b-1 ^ (a+b) 2 J n ~ 2 ^ 2 ^ 06 ( 0 + 6 - 1 ) J 


Table 7. Sequences of polynomials orthogonal over of the four major laws. 


Each measure also has an infinite sequence of monic polynomials q n (x) which are 
orthogonal with respect to that measure. They can be written as sums of Chebyshev 
polynomials of the second kind, U n (x), which satisfy U- 1 = 0, Uq{x) = 1, and 
U n (x) = 2xU n -i(x) — U n - 2 (x) for n > 1, [T2j - See Table [?] For n = 0, qo(x) = 1, 
and in general for n > 1, 

q n {x) = /3” _1 (:r- a 0 )E/ n _i ((x - oi)/(2/?i)) - /3l^~ 2 U n _ 2 ((x - a 1 )/(2j3 1 )). 

In the Wigner semicircle case the polynomials can be combined using the recursion 
rule for Chebyshev polynomials. 

4. The Wachter Law Moment Pyramid. 

Using Mathematica we can extract an interesting number pyramid from the 
Wachter moments, see Figure [T] Each triangle in the pyramid is formed by tak¬ 
ing the coefficients of a and b in the «-th Wachter moment, with the row number 
within the pyramid determined by the degree of the corresponding monomial in 
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a and b. All factors of (a + b) are removed from the numerator and denominator 
beforehand and alternating signs are ignored. 

Furtheremore, there are many patterns within the pyramid. The top row of each 
triangle is a list of Narayana numbers, which sum to Catalan numbers. The bottom 
entries of each pyramid are triangular numbers. The second-to-bottom entry on the 
right of every pyramid is a sum of consecutive triangular numbers. The second to 
both the left and right on the top row of every triangle are also triangular numbers. 


5. Moments Build Nearly-Toeplitz Jacobi Matrices. 

This section is concerned with recovering a probability distribution from its Jacobi 
parameters, cq and /3* such that they are “Nearly Toeplitz,” i.e. there exists a k 
such that for i > k all on are equal and all /3j are equal. Note that i ranges from k to 
oo. The Jacobi parameters are found from a distribution by the Lanczos iteration. 

We now state the continuous Lanczos iteration, replacing the matrix A by the 
variable x and using p = pws, Pmp, Pm, Pw to compute dot products. A good 
source is m- For a given measure p on an interval /, let 


(p(x),q(x)) = J p(x)q(x)dp, 

and ||p(a;)|| = \J(p(x),p(x)). Then the Lanczos iteration is described by Table [ij 




Lanczos on Measure p 

0 -1=0 

q-i(x) = 

= 0, q 0 (x) = 1 

for n = 

0 ,1,2 ,... 

do 

v(x) 

= xq n (x) 


— 

{q n (x),v 

(x)) 

v(x) 

= v(x) — 

fin-iqn-xix) - a n q n (x) 

Pn = 

IK*) II 


Qn +1 

(x) = v(x)/Pn 

end for 




Table 8. The Lanczos iteration produces the Jacobi parameters in a and 0. 


There are two ways to compute the integrals numerically. The first is to sample 
x and q n {x) at many points on the interval of support for qo(x) = 1 and discretize 
the integrals on that grid. The second can be done if you know the moments of p. 
If r(x) and s(x) are polynomials, (r(x),s(x)) can be computed given p's moments. 
Since the q n {x) are polynomials, every integral in the Lanczos iteration can be done 
in this way. In that case, the q n (x) are stored by their coefficients of powers of x 
instead of on a grid. Once we have reached k iterations, we have fully constructed 
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1 



1 

3 

1 

1 


3 

1 


6 

6 

4 


6 


6 

10 


20 

5 


20 

10 

10 

10 


20 

10 


60 


45 

50 


75 

15 


50 

15 

20 


1 15 

50 

50 

15 1 

21 

140 

210 

84 7 

105 

315 

189 

21 


175 

210 

35 


105 

35 




21 



21 

105 

175 

105 

21 


280 

700 

560 

140 

196 

980 

1176 

392 

28 


490 

1176 

588 

56 


490 

490 

70 




196 

56 




28 




1 28 196 490 490 196 28 1 

36 504 1890 2520 1260 216 9 

336 2520 5040 3360 720 36 

1176 4704 4704 1344 84 

1764 3528 1512 126 

1176 1008 126 

336 84 

36 

1 36 336 1176 1764 1176 336 36 1 

45 840 4410 8820 7350 2520 315 10 

540 5670 17010 18900 8100 1215 45 

2520 15120 25200 14400 2700 120 

5292 17 640 15 120 3780 210 

5292 9072 3402 252 

2520 1890 210 

540 120 

45 


Fig. 1. A number pyramid from the coefficients of the Wachter law moments. 
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Algorithm: Compute Measure from Nearly Toeplitz Jacobi Matrix. 


(1) Nearly Jacobi Toeplitz Representation: Run the continuous Lanczos algorithm 
up to step k, after which all cq are equal and all 8i are equal, or very nearly so. 
If they are equal, this algorithm will recover d/i exactly, otherwise it will find 
it approximately. The Lanczos algorithm may be run using a discretization of 
the measure /x, or its initial moments. 


(o:o:oo, Po-.oo) = Lanczos {dfi(x)). 


(2) Cauchy transform: evaluate the finite continued fraction below on the interval 
of x where it is imaginary. 


g{x) = 


1 


x — ao 


x — ai — 


pi 


Pl -2 


Oik-1 — 


m-i 


' - oik + a/( a k - x) 2 - 4:81 


(3) Inverse Cauchy Transform: divide the imaginary part by —7r, to compute the 
desired measure. 

dn(x) = --Im(g(x)). 

7r 


Table 9. Algorithm recovering or approximating an analytic measure by Toeplitz matrices with 
boundary. 


the infinite Jacobi matrix using the first batch of /x’s moments, or a discretization 
of /x. 

Step 1 can start with a general measure in which case Step 3 finds an approx¬ 
imate measure with a nearly Toeplitz representation. Step 1 could also start with 
a sequence of moments. It should be noted that the standard way to go from mo¬ 
ments to Lanczos coefficients uses a Hankel matrix of moments and its Cholesky 
factorization ([5], (4.3)). 

As an example, we apply the algorithm to the histogram of the eigenvalues of 
(X/y/m + /iiy(X/\Jrn + /x/), where X is to x n, which has Jacobi parameters cq 
and 8i that converge asymptotically and quickly. We smooth the histogram using a 
Gaussian kernel and then compute its Jacobi parameters. The reconstruction of the 
histogram is in Figure [3] We also use the above algorithm to reconstruct a normal 
distribution from its first sixty moments, see Figure [4] 

The following theorem concerning continued fractions allows one to stably re¬ 
cover a distribution from its Lanczos coefficients ai and 8i- As we have said, if the 
first batch of /r’s moments are known, we can find all ai and 8i from i = 0 to oo 


11 
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using the continuous Lanczos iteration. 


Theorem 5.1. Let p be a measure on interval I Cl with Lanczos coefficients at 
and fli, with the property that all on are equal for i > k and all are equal for 
i > k. We can recover I = \a k — 2j3 k ,a k + and we can recover dffix) using a 
continued fraction. This theorem combines Theorems 1.97 and 1.102 of m- 


9(x) 


x — ao — 


x — OL\ — 


Pi 


pi 


2/?Li 


- a k + y/{a k - x) 2 


4/3 


2 

k 


dffix) = --Im( 5 (x)). 

7r 

Figure [2] illustrates curves recovered from random terminating continued frac¬ 
tions g(x) such that the /?* are positive and greater in magnitude than the a,. In 
both cases, the above theorem allows correct recovery of the at and (which is 
not always numerically possible). In the first one, k = 5, in the second, k = 3. 

If X is an m x n, m < n matrix of normals for m and n very large, ( X/y/m + 
pi)*(X/ffirn+fj.1) has and /3, which converge to a constant, making its eigenvalue 
distribution recoverable up to a very small approximation. See Figure [3] 

We also tried to reconstruct the normal distribution, whose Jacobi parameteri¬ 
zation is not at all Toeplitz, and which is not compactly supported. Figure [4] plots 
the approximations using 10 and 20 moments. 


6. Direct computation of the Wachter law moments. 

While the moments of the Wachter law may be obtained in a number of ways, 
including expanding the Cauchy Transform, or applying the inobius inverse formula 
to the free cumulants, in this section we show that a direct computation of the 
integral is possible. 


Theorem 6.1. We find the moments of the Wachter law, m k . 

2j+4 


m k = 


a + b 


k -2 

- (« + b) 

3=0 


y/a{a + b — 1)\ 

. / 


N. 


3+1 


a(a + 6—1) 


Proof. We start by integrating the following expression by comparing it to the 
Marclienko-Pastur law. 

I rv+ _ 

Ji = — J x k yj(P+ - x) (x - p-)dx. 
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Fig. 2. Recovery from of a distribution from random oli and ft using Theorem 5.1. On top we use 
k = 5, on bottom we use k = 3. 


If x = su, dx = sdu and this integral becomes 



To compare this expression to the Marchenko-Pastur law, we need to pick s and A 
such that ^ = (1 + VX) 2 and ^ = (1 — \/A) 2 for A > 1. There are more than one 


13 
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Fig. 3. Eigenvalues taken from (X/y/m fil) 1 (X/y/m -\- fil), where X is m x n, m = 10 4 , n = 3m, 
/x = 5. The blue bar histogram is taken using hist.m, a better one was taken by convolving 
the data with Gaussian kernel. That convolution histogram was used to initialize the continuous 
Lanczos algorithm which produced five a’s and /5’s. They were put into a continued fraction as 
described above, assuming a.i and to be constant after i = 5. The continued fraction recreated 
the histogram, which is the thick red line. 


choices of each parameter, but we pick y s = ^ ( v //2+ — y/t 1 -) an d 


/^ = yy+ + yy^ 

yjy-yj^' 

Using the Narayana numbers, and the formula for the moments of the Marchenko- 
Pastur law, the integral equals 


Ji = (^{yjy-yjn)^ + iv fc+ i ^ 


ym + yjn 
yjy-yj^ 


Using a and 6, this becomes 

y a(ci + 6 - 1 ) 


Ji — 


a + b 


N, 


k +1 


a(a + b — 1 ) J 


2k+i 
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Fig. 4. The normal distribution’s Jacobi matrix is not well approximated by Toeplitz plus bound¬ 
ary, but with sufficiently many moments good approximations are possible. The above graph shows 
the normal distribution recovered by the method in this paper using 10 and 20 moments. The 
thick line is the normal computed by e~ x / 2 \/27r, and the thin lines on top of it use our algorithm. 


We also need to integrate 

T 1 r + vv+ - x )( x - m-) 

Let su = x — 1. sdu = dx. This becomes 


Ji = - 


s 

2-7T 



which by symmetry is 


s 

f * v 


11 


2 nj 


u 



Using the same technique as previously, y/s = | (y/1 — — sj\ — fi + ) and 

~~ M+ + -\/l — H- 

V 1 - V- ~ V 1 - M+ ’ 

Using the fact that the Marchenko-Pastur law is normalized, the answer is 

h = s = i (v 71 - 
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Now we are ready to find the moments of the Wachter law. Using the geometric 
series formula, 


m k 


a + b f M+ x k V(m + - x){x - H-) dx 
2t r 1 - x 

j x 3 yj{ii + - x)(x - n_)dx 


a + b 

2tt 

a 

a + b 


r + V(/d+ -x)(x~n-) 

4 - 1 _ x 


k—2 

(a + b)Y^ 

3=0 


yja(a + 6-1) 


a + b 



x 3 \/(A < + - x ){ x 

_ 6 _V. 

a(a + b — 1) J 


H-)dx 


□ 
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